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In  this  paper,  a  novel  approach  for  battery  state  of  charge  (SOC)  estimation  in  electric  vehicles  (EVs)  based 
on  an  adaptive  switching  gain  sliding  mode  observer  (ASGSMO)  has  been  presented.  To  design  the 
ASGSMO  for  the  SOC  estimation,  the  state  equations  based  on  a  battery  equivalent  circuit  model  (BECM) 
are  derived  to  represent  dynamic  behaviours  of  a  battery.  Comparing  with  a  conventional  sliding  mode 
observer,  the  ASGSMO  has  a  capability  of  minimising  chattering  levels  in  the  SOC  estimation  by  using  the 
self-adjusted  switching  gain  while  maintaining  the  characteristics  of  being  able  to  compensate  modelling 
errors  caused  by  the  parameter  variations  of  the  BECM.  Lyapunov  stability  theory  is  adopted  to  prove  the 
error  convergence  of  the  ASGSMO  for  the  SOC  estimation.  The  lithium-polymer  battery  (LiPB)  is  utilised 
to  conduct  experiments  for  determining  the  parameters  of  the  BECM  and  verifying  the  effectiveness  of 
the  proposed  ASGSMO  in  various  discharge  current  profiles  including  EV  driving  conditions  in  both  city 
and  suburban. 

©  2013  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

Due  to  the  progressive  rise  of  petrol  costs  and  environmental 
concerns  on  exhaust  emission  from  petrol-driven  vehicles,  the  eco- 
friendly  electric  vehicles  (EVs)  have  greatly  exhibited  promising 
potential  to  revive  as  the  leading  means  of  transportation  in  up¬ 
coming  decades.  The  EV  performance  is  highly  dependent  on  bat¬ 
tery  characteristics  such  as  operation  voltage,  temperature,  charge 
or  discharge  rate  and  aging.  Among  existing  battery  types  appli¬ 
cable  to  EVs,  lithium-polymer  battery  (LiPB)  is  widely  recognised  as 
the  most  capable  candidate  for  the  development  and  innovation  of 
the  new  generation  EVs.  Compared  with  other  types  of  EV  batteries 
such  as  lead-acid  batteries,  nickel-cadmium  batteries  and  nickel— 
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metal  hydride  batteries,  LiPBs  are  superior  in  terms  of  high  energy 
and  power  density,  broad  operating  temperature  range,  rapid 
charge  capability,  no  memory  effects,  long  cycle  life  and  extremely 
low  self-discharge  rate  [1,2]. 

The  amount  of  battery  available  capacity  is  closely  related  to  the 
state  of  charge  (SOC),  which  is  considered  as  one  of  the  key  factors 
in  battery  management  system  (BMS)  for  supporting  optimal  bat¬ 
tery  performance  and  safety  in  EVs.  The  accurate  battery  SOC 
indication  is  essential  for  predicting  a  reliable  travelling  range, 
maximising  the  efficiency  of  battery  energy  utilisation  and  pre¬ 
venting  the  batteries  in  EVs  from  over-charging  or  over¬ 
discharging.  Unfortunately,  the  SOC  involves  in  intrinsic  electro¬ 
chemical  processes  of  a  battery,  and  it  cannot  be  directly  measured 
by  a  sensor.  It  should  be  estimated  by  an  advanced  mathematical 
algorithm  with  the  aids  of  measurable  signals  such  as  terminal 
voltage  and  current  from  the  battery. 
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A  number  of  techniques  have  been  proposed  to  estimate  the  SOC 
of  a  battery  and  each  one  has  its  own  advantages  in  certain  aspects 
[3-5].  The  most  straightforward  approach  is  the  ampere-hour  (Ah) 
counting,  which  simply  integrates  the  battery  charge  and  discharge 
currents  over  time.  The  SOC  can  be  calculated  by  referring  to  the 
calibration  point  at  the  fully  charged  battery.  This  method  is  inex¬ 
pensive  to  be  implemented  in  hardware,  but  the  estimation  accu¬ 
racy  is  strongly  dependent  on  the  sampling  frequency  and  the 
precision  of  the  current  sensor.  Moreover,  this  non-model  and 
open-loop-based  estimator  can  easily  accumulate  errors  caused  by 
embedded  noises  and  current  measurement  drift  and  it  is  also 
incapable  of  determining  the  initial  SOC.  An  enhanced  version  of  Ah 
counting  has  shown  improved  SOC  estimation  results  by  online 
evaluating  Coulombic  efficiency  with  recalibration  of  the  battery 
capacity  [6]. 

The  impedance  measurement  has  been  declared  as  an  effective 
technique  for  the  SOC  estimation  [7—9].  It  loads  a  series  of  small 
amplitude  a.c.  signals  into  a  battery  to  detect  the  responses  of  the 
battery  in  the  wide  range  of  frequencies.  The  SOC  can  be  obtained 
by  analysing  the  battery  impedance,  but  the  measured  impedances 
cannot  completely  represent  the  dynamics  of  batteries  in  the  case 
of  large  discharge  current  in  EVs.  Another  drawback  of  this 
approach  is  that  a  set  of  the  bulky  and  costly  auxiliary  equipment 
such  as  the  signal  generator  and  the  impedance  spectroscope  is 
required  to  carry  out  impedance  measurement. 

Another  category  of  the  SOC  estimation  methods  is  relied  on 
computational  and  intelligence-based  strategies,  which  encompass 
artificial  neural  networks  (NNs),  fuzzy  neural  networks,  adaptive 
fuzzy  neural  networks  and  support  vector  machine  [10-13].  These 
data-oriented  approaches  can  accurately  estimate  the  SOC  for  all 
kinds  of  batteries  in  the  absence  of  the  details  of  batteries,  but  they 
require  a  large  number  of  training  sample  data  to  train  the  NNs. 
Therefore,  they  demand  more  powerful  and  costly  data  processing 
chips  to  handle  the  massive  computation  loads  in  the  BMSs. 
Furthermore,  the  SOC  would  be  unpredictable  in  case  of  discharge 
current  profiles  loaded  in  EV  batteries  deviated  from  those  repre¬ 
sented  by  the  training  data. 

The  Kalman  filter  (KF),  as  a  classical  state  estimation  method  for 
dynamic  systems,  was  also  developed  to  estimate  the  SOC  based  on 
a  linear  model  [14,15].  For  the  battery  represented  by  a  nonlinear 
model,  some  advanced  KF  techniques  such  as  extended  KF  (EKF) 
[16,17],  sigma-point  KF  (SPKF)  [18,19]  and  unscented  KF  (UKF) 
[20,21]  were  proposed  to  achieve  online  SOC  estimation.  The 
essential  idea  of  the  EKF  approach  is  to  transform  a  nonlinear 
system  into  a  linear  system  by  linearising  the  nonlinear  function 
based  on  the  first  order  Taylor  series  expansion,  such  a  process 
gives  rise  to  large  linearisation  errors  and  complicated  computation 
of  the  Jacobian  matrix  which  may  lead  to  the  instability  of  the  filter 
and  inaccurate  estimation  for  highly  nonlinear  battery  systems  in 
EVs.  Instead  of  the  local  linearisation  in  the  EKF,  the  SPKF  and  UKF 
approaches  use  an  unscented  transformation  to  approximate  a 
Gaussian  distribution  of  the  state  random  variable  with  a  set  of 
sample  points  or  sigma  points  and  offer  better  SOC  estimation  re¬ 
sults  in  terms  of  accuracy  and  robustness  [18-21]. 

All  these  KF-based  SOC  estimation  algorithms,  however,  require 
accurate  battery  model  parameters  with  the  assumption  that 
constant  values  of  the  process  and  measurement  noise  covariance 
are  priori  known,  which  are  hardly  practical  and  error-prone  due  to 
the  complex  electrochemical  reactions  inside  batteries  for  EV 
driving  conditions.  Furthermore,  the  constant  values  of  noise 
covariance  can  result  in  remarkable  errors  and  divergence  in  the 
battery  SOC  estimation.  Later,  both  an  adaptive  EKF  (AEKF)  and  an 
adaptive  UKF  (AUKF)  have  been  developed  to  estimate  the  values  of 
the  process  and  measurement  noise  covariance  during  the  online 
process  [22,23].  They  have  demonstrated  better  precision  in  the 


SOC  estimation  and  filter  divergence  restraint  at  the  expense  of 
higher  complexity  and  computational  cost. 

More  recently,  the  Hx  observer-based  method  has  been  pro¬ 
posed  to  estimate  the  SOC  without  the  requirement  of  the  exact 
statistical  properties  of  the  battery.  This  method  minimises  the 
errors  of  system  and  measurement  so  that  the  SOC  estimation  error 
is  less  than  a  given  attenuation  level  [24,25],  where  an  alternative 
feedback  gain  is  employed  to  tackle  modelling  errors  and  distur¬ 
bances.  Similarly,  the  sliding  mode  observer  (SMO)  based  SOC 
estimation  method  has  been  adopted  to  overcome  the  un¬ 
certainties  of  battery  model,  external  disturbances  and  measure¬ 
ment  noises  [26,27].  Nevertheless,  this  method  relies  on  the 
exhaustive  understanding  of  battery  dynamics  for  the  appropriate 
selection  of  the  SMO  parameters  such  as  uncertainty  boundaries 
and  switching  gains,  leading  to  the  trade-off  between  the  chatter¬ 
ing  magnitude  and  the  convergence  speed  in  the  SOC  estimation. 

In  this  paper,  a  novel  approach  for  the  SOC  estimation  based  on 
an  adaptive  switching  gain  sliding  mode  observer  (ASGSMO)  has 
been  proposed.  Comparing  with  constant  switching  gains  SMO,  the 
ASGSMO  is  able  to  dynamically  adjust  the  switching  gains  in 
response  to  the  tracking  errors,  and  guarantee  the  reachability  of 
sliding  mode  surface  and  trigger  the  sliding  mode.  Once  the  sliding 
mode  is  activated,  the  switching  gains  are  self-tuned  to  “adequate” 
levels  to  counteract  the  modelling  errors  and  reduce  the  chattering 
magnitudes,  thereby  improving  the  SOC  estimation  accuracy. 

The  remaining  part  of  this  paper  is  organised  as  follows.  In 
Section  2,  a  battery  equivalent  circuit  model  (BECM)  is  presented  to 
characterise  the  discharge  behaviours  of  the  LiPB  in  the  presence  of 
parameter  uncertainties.  The  detailed  procedures  to  identify  the 
BECM  parameters  are  also  explained  in  this  section.  In  Section  3, 
the  ASGSMO  design  methodology  is  elaborated  for  the  SOC  esti¬ 
mation.  The  proposed  ASGSMO  is  validated  for  SOC  estimation  by 
experimental  results  in  Section  4,  followed  by  conclusions  in  Sec¬ 
tion  5. 

2.  Battery  modelling 

2.1.  Battery  equivalent  circuit  model 

There  have  been  numerous  attempts  to  establish  a  precise  bat¬ 
tery  model  for  achieving  an  accurate  battery  state  estimation.  The 
battery  equivalent  circuit  models  consisting  of  circuit  components 
such  as  capacitors,  resistors,  diodes  and  voltage  sources  have  been 
widely  studied  and  developed  to  capture  dynamic  characteristics  of 
a  battery  for  reducing  modelling  errors  [28-30].  These  circuit- 
based  models  are  also  applied  to  the  battery  SOC  estimation  due 
to  their  state  equations  are  intuitively  derived  from  circuit  analysis 
for  mathematical  computation.  In  this  paper,  without  the 
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Fig.  1.  Schematic  diagram  of  BECM. 
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consideration  of  a  hysteresis  effect,  a  BECM  is  used  to  model  the 
LiPB  dynamics  and  its  schematic  diagram  is  shown  in  Fig.  1,  where 
the  variations  of  the  circuit  parameters  are  also  included  to 
represent  the  modelling  errors. 

In  the  BECM,  a  capacitor  Cn  with  the  nominal  capacity  represents 
the  total  stored  energy  in  the  battery.  Vt  and  I  represent  the  battery 
terminal  voltage  and  discharge  current,  respectively.  The  voltage 
across  Cn  as  Vsoc  varies  in  the  range  of  the  SOC  from  0%  (the  fully 
discharged)  to  100%  (the  fully  charged),  and  it  represents  the  SOC 
quantitatively.  A  voltage-controlled  current  source  /  (Vsoc)  and 
voltage-controlled  voltage  source  Voc  (Vsoc)  are  used  to  bridge  a 
nonlinear  relationship  between  the  SOC  (i.e.,  Vsoc)  and  the  open 
circuit  voltage  (OCV)  of  the  LiPB  as  shown  in  Fig.  2,  which  can  be 
derived  by  fitting  the  experimental  data  from  a  pulse  current 
discharge  (PCD)  profile  (see  details  later). 

Two  parallel  resistive  and  capacitive  branches  consist  of  the 
electrochemical  polarisation  resistance  Rpe,  capacitance  Cpe,  the 
concentration  polarisation  resistance  Rpc  and  capacitance  Cpc. 
which  reflect  the  short-term  and  long-term  transient  responses 
due  to  the  relaxation  effect  of  battery,  respectively.  This  relaxation 
effect  represents  the  slow  convergence  of  battery  terminal  voltage 
to  the  OCV  (Voc)  at  its  equilibrium  state  after  several  hours  of 
relaxation  at  the  end  of  charging  or  discharging.  It  is  mostly  caused 
by  a  diffusion  effect  and  a  double-layer  charging  or  discharging 
effect  in  the  battery  [28].  An  ohmic  resistor  characterises  the 
charging  or  discharging  energy  losses  of  battery.  The  self-discharge 
resistance  is  neglected  in  this  model  as  the  LiPB  has  extremely  low 
self-discharge  rate. 

The  symbols  AV0C,  A ftin,  Aftpe,  ACpe,  A Rpc  and  ACpc  associated 
with  corresponding  circuit  components  represent  variations  of  the 
circuit  parameters  caused  by  the  combination  of  modelling  errors 
and  disturbances.  There  are  two  main  causes  for  modelling  errors. 
Firstly,  the  circuit  parameters  of  the  BECM  are  determined  as  the 
constant  values,  but  in  fact  they  are  varying  with  the  battery  SOC 
and  temperature  in  different  discharge  current  profiles.  Secondly, 
the  nonlinear  relationship  of  OCV— SOC  curve  is  piecewise  line¬ 
arised  as  derived  later  and  it  results  in  large  linearisation  errors. 

The  SOC  is  a  value  between  0%  and  100%  that  indicates  the  ratio 
of  available  capacity  to  the  nominal  capacity  in  the  battery.  It  can  be 
obtained  through  the  Ah  counting  approach  by  taking  a  time  in¬ 
tegral  of  measured  battery  current  as  expressed  by  Eq.  (1). 


u 

S0C(f,)  =  Soc(to)  -  (1  /Cn)  J  Jj/(T)dT  (1) 

to 

where  Soc(to)  denotes  the  initial  SOC,  Cn  is  the  battery  nominal 
capacity  in  ampere-hours,  /( t)  is  the  instantaneous  discharge  cur¬ 
rent  (assumed  to  be  positive),  and  p  is  the  Coulombic  efficiency, 
which  is  usually  taken  as  1  for  discharging  and  smaller  than  or  close 
to  1  for  charging  for  the  LiPBs  in  the  wide  range  of  current  and 
temperature. 

According  to  the  Kirchhoff  voltage  law  in  circuit  theory,  the 


battery  terminal  voltage  Vt  can  be  determined  by  Eq.  (2). 

Vt  —  V0C(S0C)  —  Vpe  —  Vp  c  —  IRin  +  A uncertain  (2) 

The  time  derivatives  of  SOC  and  polarisation  voltages  give 

Soc  =  -(//C„)+  A/2  (3) 

Vpe  =  —Vpe/(f^peCpe)  + //Cpe  +  A/3  (4) 

Vpc  =  —  Vpe/  (RpcCpc)  +  //Cpc  +  A/4  (5) 


where  Vpe  and  Vpc  denote  the  electrochemical  and  concentration 
polarisation  voltages  across  Cpe  and  Cpc,  respectively.  The  terms 
A uncertain ,  A f2,  A/3  and  A/4  not  only  represent  model  uncertainties 
caused  by  modelling  errors,  but  also  time-varying  internal  and 
external  disturbances. 

Despite  the  nonlinearity  of  OCV— SOC  curve  as  shown  in  Fig.  2, 
there  exists  a  piecewise  linear  relationship  between  the  OCV  and 
the  SOC  in  a  certain  range  of  SOC  indicated  by  the  red  dots  in  the 
curve.  Therefore,  the  OCV  can  be  expressed  as  a  function  of  the  SOC 
by  using  piecewise  linearisation  method. 

Voc  (Soc)  =  kSoc-\-v  (6) 

where  the  values  of  k  and  v  are  varying  in  different  SOC  ranges. 

Due  to  the  fast  sampling  rate  [27],  the  changing  rate  of  discharge 
current  I  can  be  negligible,  namely  d//dt  =  0,  the  time  derivative  of 
terminal  voltage  Vt  in  Eq.  (2)  with  the  substitutions  of  Eqs.  (3)— (6) 
gives 

Vt  =  — K(I/Cn)  +  Vpe/ (RpeCpe)  ~  //Cpe  +  Vpc/ (RpcCpc)  — //Cpc  +  A/| 

(7) 

Solving  /  in  Eq.  (2)  and  substituting  it  into  Eq.  (3)  as  well  as  re¬ 
arranging  Eqs.  (3)— (6)  result  in  the  state-space  equations  of  the 
BECM  as  follows 

Vt  =  Vt  +  0\  Vqc(Soc)  —  G3 Vpe  —  Q^Vpc  —  b\l  +  A/l 

SOC  =  CI2  Vt  —  Cl2V0c(S0C)  +  C^Vpe  +  02Vpc  +  A/2  tgA 

Vpe  =  —0/[Vpe  +  b2 /  +  A/3 

Vpc  =  —  O3  VpC  +  b^I  +  A/4 

where  0\  =  1  l(RpeCpe)  +  V(^pcCpc)*  ^2  —  l/(^inQi),  ^3  =  'Ll  (RpcCpc), 
O4  —  'Ll  (RpeCpe),  hi  =  K/Cn  +  Rinl(RpeCpe)  +  1  /Cpe  +  Rinl(RpcCpc)  +  V^pc* 
/>2  —  'Ll  Cpe  and  =  J/Cpc. 

If  the  input  and  the  output  of  circuit  model  are  defined  as  u(t)  =  / 
and y(t)  =  Vt,  respectively,  and  the  system  state  variables  are  chosen 
as  Vt,  Vqc(Soc),  Vpe  and  Vpc,  then  Eq.  (8)  of  the  BECM  can  be  concisely 
expressed  in  matrices 

x(t)  =  Ax(t)  +  Bu(t )  +  A/(x,  u,  t) 


(9) 
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y(t)  =  Cvm  (10) 

where 


"-«i 

-a3 

-a4- 

Mr 

A  = 

-02 

Cl2 

Cl2 

,  B  = 

0 

0 

0 

-a4 

0 

b2 

0 

0 

0 

~a3- 

-  b3  _ 

C=[l  0  0  0]  andx(t)  =  [VtV0c(Soc)VpeVpc\T . 

The  unknown  function  A /(x,u,t)  represents  the  matched  model 
uncertainties  and  it  can  be  assumed  that 

A f{x,u,t)  =  r%{x,t)  (ii) 

where  r  is  the  model  uncertainties  input  matrix  and  the  function 
?(x,t)  is  unknown  but  bounded 

|?(x,t)|<^  VxsR4,  t  >  0  (12) 

where  the  bound  ^  is  a  positive  constant  and  can  be  determined  by 
the  largest  modelling  errors  between  the  experimental  data  and 
the  BECM  of  the  LiPB.  According  to  Eq.  (11),  Eq.  (9)  is  rewritten  as 

x(t)  =  A*(t)+Bu(t)  +  r?(x,t)  (13) 

The  success  of  applying  the  ASGSMO  to  estimate  the  internal 
states  of  the  battery  depends  on  the  stability  and  observability  of 
the  corresponding  linear  parts  of  the  BECM. 

Remark  1:.  The  eigenvalues  of  the  BECM  system  are  obtained  by 


det  (AJ0  -A)  =  det 


A  +  Oq 

fli 

a3 

-a4 

t*2 

A  +  a2 

&2 

«2 

0 

0 

A  +  a4 

0 

0 

0 

0 

A  +  a3 

=  0 


(14) 


where  Jo  is  an  identity  matrix  with  appropriate  dimensions.  The 
eigenvalues  of  the  matrix  A  are  X\  =  -1  /(ftpeCpe)  -  1  /(RpcCpc),  A2  =  0, 
A3  =  -l/(RpeCpe)  and  A4  =  -l/(RpcCpc).  Therefore,  based  on  the 
modern  control  theory,  the  BECM  system  is  stable. 

Remark  2:.  The  observability  matrix  of  the  BECM  system  is 

1  0  0  0  " 

— Q\  Cl]  Q3  Cl/\ 

Oft (3,  1)  Oft (3,  2)  Oft(3,  3)  Oft(3,  4) 

Oft (4,  1)  Oft (4,  2)  Oft (4,  3)  0„(4,4)_ 

(15) 

where 


Ob(A,C)  = 


C 
CA 
CA 2 
G43 


Oft (3,  1)  =  a\  +  a2di,  Oft (3,  2)  =  -a\-a2a-l, 


0&(3,  3)  =  a3ai  +  a2ai  +  a3a4, 
Ob( 3,  4)  =  a4ci]  +a2a!  +  a3a4, 

Oft (4,  1)  =  -a3  -  2a2a3  -  a2a3, 
Oft  (4,  2)  =  a3  +  2a2a3  +  a2a\, 


Oft (4, 3)  =  -a3a^-a2a3a1  -a2aj-a^aj -a4a3a1-a2a4a-i  -a|a3, 


Oft(4, 4)  =  -a4a3-a2a4ai  - a2a\-a2a i  -a4a3ai  -a2a3a i  -a3a4. 


On  the  conditions  of  l/(ftpeCpe)  =£  1,  l/(ftpcCpc)  =£  1  and  1/ 
( RinCn )  =£  1,  which  are  always  satisfied  for  the  battery  in  normal 
operating  conditions,  then  it  enables  the  full  rank  of  Eq.  (15)  so  that 
the  BECM  system  is  observable  and  thus  the  internal  states  of  the 
battery  can  be  estimated  by  an  observer. 

The  advantages  of  this  BECM  are  that  it  relatively  captures  the 
dynamic  behaviours  of  a  LiPB  with  reasonable  computational  cost. 
Additionally,  the  modelling  errors  and  disturbances  also  can  be 
compensated  by  the  robust  ASGSMO. 


2.2.  Battery  model  parameters  determination 

The  circuit  parameters  of  the  BECM  in  Fig.  1  are  obtained  from 
the  transient  response  of  terminal  voltage  by  executing  the  PCD 
profile  as  a  battery  characterisation  test  at  room  temperature  [27]. 
The  Turnigy  LiPB  in  Fig.  3  is  used  in  the  experiments.  It  comprises  of 
a  LiMn204  cathode  and  an  artificial  graphite  anode,  which  are 
usually  laminated  with  the  solid  electrolyte  and  the  separator 
sheets  and  enclosed  in  a  foil-type  case.  The  solid  electrolyte  enables 
batteries  safer  and  allows  the  use  of  higher  energy  electrode  ma¬ 
terials  designed  for  high  power  applications.  The  LiPB  has  a  nom¬ 
inal  capacity  of  5.0  Ah  and  a  nominal  voltage  of  3.7  V.  The 
dimension  of  the  cell  is  135  mm  x  50  mm  x  9  mm  and  the  weight 
of  the  cell  is  130  g.  A  battery  testing  platform  as  illustrated  in  Fig.  4 
is  developed  to  conduct  experiments,  which  consists  of  a  host 
computer,  a  programmable  power  supply  (Sorensen  DLM50-60),  an 
electronic  load  (Prodigit  3320)  and  a  safety  control  switch  box.  The 
battery  testing  platform  can  control  charge/discharge  battery  cur¬ 
rent,  sample  battery  terminal  voltage  current  and  temperature,  and 
then  record  the  data  into  the  computer  via  a  graphic  user  interface 
program  designed  by  using  the  LABVIEW  software. 

The  PCD  profile  as  an  off-line  method  was  performed  for  the 
BECM  parameters  identification  as  shown  in  Fig.  5.  Initially,  the  LiPB 
was  fully  charged  (SOC  =  100%)  at  room  temperature,  and  left  it  in 
the  open  circuit  state  for  an  hour  so  that  the  OCV  is  able  to  reach  its 
equilibrium  state.  Then,  the  LiPB  was  discharged  by  a  sequence  of 
5.0  A  pulse  currents  consisting  of  6-min  discharge  and  1  h  rest 
between  two  consecutive  pulse  currents  to  allow  the  OCV  return  to 
the  equilibrium  state.  Each  pulse  current  discharges  approximately 
10%  of  the  nominal  capacity,  which  is  equivalent  to  10%  of  the  SOC 
reduction.  The  procedure  of  the  pulse  current  discharge  and  re¬ 
covery  was  iterated  until  the  terminal  voltage  reached  the  cut-off 
voltage  of  2.7  V  as  the  battery  fully  discharged  state  (SOC  =  0%). 
The  corresponding  voltage  response  is  shown  in  Fig.  5.  It  can  be 
seen  that  totally  ten  sets  of  transient  response  in  terminal  voltage 
generated  to  determine  circuit  parameters  for  ten  different  SOCs. 
Meanwhile,  with  the  aid  of  Ah  counting  method,  the  nonlinear 
relationship  between  the  OCV  and  the  SOC  over  entire  operation 
range  can  be  obtained  by  conducting  the  PCD  profile  on  the  LiPB 
and  the  corresponding  nonlinear  curve  as  shown  in  Fig.  2. 

Intuitively,  during  relaxation  period  each  terminal  voltage  can 
be  represented  by 

Vt(t)  =  Voc  —  Vpeexp(  —  t/xpe)  —  VpCexp(  —  t/ipc)  (16) 


Fig.  3.  Photo  of  Turnigy  LiPB. 
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Fig.  4.  Testing  platform  for  the  LiPB. 


where  ipe  and  ipc  are  the  time  constants  for  the  electrochemical 
polarisation  and  concentration  polarisation  voltages  during  tran¬ 
sient  response,  respectively.  Voc  is  the  OCV  after  a  full  relaxation.  Vpe 
and  Vpc  are  the  voltages  of  the  electrochemical  and  concentration 
polarisation  capacitors,  respectively.  According  to  the  circuit  anal¬ 
ysis,  the  following  equations  are  derived  to  calculate  the 


parameters 

Rpe  —  Vpe/J 

(17) 

Rpc  =  Vpc/I 

(18) 

Cpe  —  ^pe/Rpe 

(19) 

Cpc  —  ^pc/Rpc  (20) 

The  ohmic  resistance  is  given  by 

Rin  =  A  Vt/I  (21) 

As  an  example,  the  transient  response  of  terminal  voltage  at  the 
second  pulse  current  indicated  by  a  red  circle  in  Fig.  5  was  used  to 
extract  the  circuit  parameters,  where  the  circled  section  is  magni¬ 
fied  in  Fig.  6.  With  the  circuit  parameters  identified  by  using  a  curve 
fitting  technique  based  on  Eq.  (16)  and  calculated  by  Eqs.  (17)— (21 ), 
the  curve  fitting  results  are  plotted  in  Fig.  6  in  the  red  dash  line. 

The  pulse  current  discharge  was  repeated  ten  times  corre¬ 
sponding  to  different  SOCs,  thus  ten  groups  of  circuit  parameters 
were  obtained  and  summarised  in  Table  1  with  the  curve  fitting 
errors  represented  by  root  mean  square  errors  (RMSEs).  The  BECM 
has  a  small  deviation  from  the  LiPB  in  the  training  profile  due  to 


Fig.  5.  Pulse  current  discharges  and  corresponding  terminal  voltages. 


Fig.  6.  Transient  response  of  the  LiPB  and  BECM  for  second  pulse  current. 
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Table  1 

Circuit  parameters  of  BECM  and  RMSEs  at  different  SOCs. 


SOC 

Rin  (mfl) 

RPe  (mfl) 

Cpe  (1<F) 

Rpc  (mfl) 

CPc  (1<F) 

RMSE 

90% 

102.1 

4.62 

4.05 

1.34 

12.26 

0.0256 

80% 

102.5 

4.96 

4.93 

2.86 

14.33 

0.0293 

70% 

102.9 

4.23 

7.05 

4.92 

16.67 

0.0268 

60% 

103.0 

3.58 

5.51 

2.17 

17.52 

0.0238 

50% 

102.9 

3.62 

4.81 

1.19 

12.98 

0.0266 

40% 

103.1 

4.12 

6.03 

1.45 

13.94 

0.0272 

30% 

103.2 

3.36 

6.62 

1.63 

16.46 

0.0285 

20% 

103.3 

3.16 

6.02 

1.52 

18.71 

0.0245 

10% 

103.3 

4.13 

6.95 

3.27 

20.54 

0.0274 

0% 

104.6 

10.31 

6.76 

12.23 

21.387 

0.0289 

non-zero  RMSEs.  Since  the  parameters  at  80%  of  SOC  caused  the 
largest  fitting  error  as  listed  in  Table  1,  they  were  used  to  calculate 
the  system  parameters  in  Eq.  (8). 

To  illustrate  the  modelling  errors  caused  by  using  constant  cir¬ 
cuit  parameters  and  demonstrate  the  robustness  of  the  ASGSMO  in 
presence  of  the  modelling  errors,  the  constant  current  discharge 
(CCD),  city  and  suburban  current  profiles  are  concurrently  loaded  to 
the  LiPB  and  the  BECM,  where  the  city  and  suburban  current  pro¬ 
files  simulate  more  realistic  EV  driving  conditions  in  both  city  and 
suburban,  respectively  [31].  These  profiles  with  the  average 
discharge  current  of  l/3Cn,  lCn  and  1.5Cn  are  adopted  to  conduct 
experiment  in  this  paper,  which  are  corresponding  to  normal  EV 
driving  conditions.  As  an  example,  three  profiles  with  the  average 
discharging  rate  lCn  and  their  experimental  results  in  terms  of 
battery  and  the  estimated  voltages  from  the  BECM  are  shown  in 
Figs.  7-9,  respectively.  It  can  be  seen  that  the  BECM  output  voltages 
are  not  able  to  track  the  LiPB  voltage  and  there  always  exists 
modelling  errors  in  each  profile  and  their  RMSEs  associated  with 
various  discharging  rates  are  listed  in  Table  2.  It  shows  that 
modelling  errors  are  increased  in  large  discharge  current  rates  and 
the  model  accuracies  are  relatively  low  in  such  conditions  due  to 
the  fixed  model  parameters  obtained  from  the  off-line  identifica¬ 
tion.  To  overcome  the  model  inaccuracy  issue,  the  ASGSMO  will  be 
designed  to  provide  the  robust  SOC  estimation  by  compensating 
the  modelling  errors  and  uncertainties  in  the  next  section. 


4.5  i 


LIPB  Vt  City  lCn 


2.5  -I - 1 - 1 - 1 - 1 — 

0  1000  2000  3000  4000 

Time  [s] 


Fig.  8.  City  lCn  profile  and  terminal  voltages  of  LiPB  and  BECM. 

3.  Design  of  adaptive  switching  gain  sliding  mode  observer 
for  SOC  estimation 

The  SMOs  with  constant  switching  gains  for  the  SOC  estimation 
have  demonstrated  the  robustness  to  modelling  errors  and  dis¬ 
turbances  with  properly  chosen  switching  gains  [26,27].  Since  the 
modelling  errors  and  the  SOC  estimation  errors  are  closely  related 
to  battery  discharge  current  profiles,  which  vary  between  1  /3Cn  and 
1.5Cn  discharge  rates  averagely  corresponding  to  normal  EV  driving 
conditions,  the  robust  and  accurate  SOC  estimation  is  unlikely  to 
be  achieved  by  using  SMOs  with  un-adjustable  switching  gains 
for  various  discharge  current  rates.  Either  underestimated  or 


Fig.  9.  Suburban  1  Cn  and  terminal  voltages  of  LiPB  and  BECM. 
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Table  2 

Modelling  errors  of  BECM  expressed  as  RMSEs  in  CCD/city/ 
suburban  profiles. 
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Current  profile 

RMSE  of  BECM 

CCD  1.5C„ 

0.164 

CCD  1C„ 

0.128 

CCD  1/3C„ 

0.096 

City  1.5Cn 

0.284 

City  1C„ 

0.238 

City  1/3  Cn 

0.193 

Suburban  1.5Cn 

0.185 

Suburban  lCn 

0.137 

Suburban  1/3  Cn 

0.105 

overestimated  switching  gains  can  lead  to  the  poor  tracking  per¬ 
formance  or  the  undesired  chattering  phenomena  in  the  SOC 
estimation  results,  respectively. 

For  instance,  the  ordinary  SMO  with  conservative  switching 
gains  is  able  to  estimate  the  SOC  with  low  magnitude  of  chattering 
in  1/3  Cn  CCD  profile  as  shown  in  Fig.  10,  but  its  tracking  perfor¬ 
mance  under  1.5Cn  CCD  profile  is  very  poor  with  average  SOC 
estimation  error  more  than  10%  as  shown  in  Fig.  11.  On  the  other 
hand,  the  tracking  performance  of  the  SMO  with  large  switching 
gains  is  significantly  improved  for  1.5Cn  CCD  profile  with  most  of 
estimation  errors  bounded  within  3%  as  shown  in  Fig.  12,  but  it 
produces  high  magnitudes  of  chattering  phenomena  for  the  SOC 
estimation  in  l/3Cn  CCD  profile  as  shown  in  Fig.  13,  and  significant 
chattering  effect  can  result  in  instability  of  observer  and  the  blurred 
SOC  estimation  results.  In  other  words,  the  SOC  estimation  error 
would  be  large  or  divergent  if  the  inappropriate  switching  gains 
were  used  in  the  highly  dynamic  current  profile. 

Therefore,  in  order  to  accurately  estimate  the  SOC  for  various 
discharge  current  profiles,  the  ASGSMO  based  on  the  BECM  is 
proposed  as  follows 


fc(t)  =  Ax(t)  +  Bu(t)  +  I<(y(t)  -  Cx(t ))  +  v  (22) 


Fig.  11.  Ordinary  SMO-based  SOC  estimation  with  conservative  switching  gains  for 
CCD  1.5C„. 


m  =  cm  (23) 

where  the  feedback  gain  matrix,  I<  and  the  adaptive  switching  gains 
function,  v  can  be  designed  so  that  the  stability  and  robustness  of 
the  ASGSMO  is  guaranteed. 


U  =  fO(ey(t))r sgn  (ey(t))  if  ey(t)=£ 0 

1  0  otherwise 

The  discontinuous  switching  term  is  defined  as 


(24) 
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Fig.  10.  Ordinary  SMO-based  SOC  estimation  with  conservative  switching  gains  for 
CCD  1/3C„. 
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Fig.  13.  Ordinary  SMO-based  SOC  estimation  with  large  switching  gains  for  CCD  l/3Cn. 


ey(t)  >  0 
ey(t)  <  0 


(25) 


where  ey(t)  =  y(t )  -  y(t)  =  Ce(t )  and  the  adaptive  switching  gain 
0(Ce(t))  are  updated  by  the  following  adaptive  law 


?(ey(t))  m  a\ey(t)\  (26) 

where  a:  is  a  positive  constant,  which  can  be  chosen  to  adjust  the 
speed  adaptation  for  switching  gains  function. 

The  state  estimation  errors  are  defined  as  e(t)  =  x(t)  -  x(t)  and 
subtracting  Eq.  (22)  from  Eq.  (13)  gives  the  dynamical  state  esti¬ 
mation  error  system  as 

e(t)  =  A0e(t)  +  T?(x,  t)-v  (27) 

where  A0  =  (A  -  I<C). 

According  to  the  Remark  2,  there  exists  a  system  feedback  gain 
matrix  I<  such  that  Ao  =A-KC  is  a  strictly  Hurwitz  matrix,  where  I< 
can  be  directly  obtained  by  using  either  the  pole  placement  method 
or  linear  quadratic  regulator  (LQR)  method. 

Based  on  an  appropriate  Lyapunov  equation  (LE),  there  exists  a 
symmetric  definite  matrix  P  as  the  solution  of  the  LE 

A0P  +  PAt0  =  -Q  (28) 

where  Q  is  a  symmetric  positive  definite  matrix  such  that  the 
structural  constraint 

C  =  rTP  (29) 

is  satisfied  by  a  Lyapunov  matrices  pair  (P,  Q)  for  Ao. 

It  can  be  seen  from  Eq.  (27)  that  if  the  adaptive  switching  gains 
function  v  is  properly  adjusted  corresponding  to  the  state  estima¬ 
tion  errors  so  that  the  asymptotically  stability  of  the  reconstruction 
error  system  is  guaranteed  if  Eq.  (29)  is  fullfied.  Thus,  a  sliding 
mode  motion  can  be  induced  on  the  system  error  states  and  the 
state  estimation  errors  is  asymptotically  converged  to  zeros  as  time 


tends  to  infinity,  a  Lyapunov  function  candidate  including  both 
state  errors  and  adaptive  error  is  considered  as  follows 

v(e(t),6(t))  =  (1/2)(eTPe  +  a~'d2^  (30) 

where  6(t)  =  6(t)  -  6d  is  the  adaptation  error  between  estimated 
switching  gain  and  desired  switching  gain.  Taking  the  time  deriv¬ 
ative  of  the  candidate  of  Lyapunov  function  V  gives 

V  =  (1/2)(eTPe  +  eTPe)  +  a~'l66 

=  (-l/2)[eT{A0P  +  PAToy\+(-i/2){[frT-vT]Pe 
+  eTP[r f  —  u]  |  +  a-1 66 

=  -(1  /2)eTQe  +  ?. rTPe  -  vTPe  +  a”1  p(t)  -  6dj  ~6 

=  -(1/2 )eTQe  +  £Ce  -  6(t) sgn  (ey)TTPe  +  or1  [?(f)  -  6dj 
x  («|ey(t)|) 

=  -(\/2)eTQe  +  aey-d(t)\ey(t)\+[d(t)-6d]\ey(t)\ 

=  -(1/2 )eTQe  +  [f  -  6d sgn  (ey)]ey 

(31) 

The  term  -(l/2)eTQe  is  always  negative  and  the  resultant 
equation  will  be  V  <  0  as 

'?  -  Odsgn  (ey)Jey  <  0  for  ey  >  0  if  0d  > 

-  Odsgn  (eyj  ey  <  0  for  ey  <  0  if  6d  >  \p  ^  ' 

As  |f(x,  t)\  <  and  there  exists  a  desired  finite  non-negative 
switching  gain  dd  such  that  dd  >  \p,  leading  to  V  <  0,  which  sat¬ 
isfies  the  second  method  of  Lyapunov  stability  theory.  Therefore, 
the  state  estimation  error  system  Eq.  (27)  asymptotically  converges 
to  zero  as  time  tends  to  infinity.  Once  the  sliding  surface  is  reached, 
the  influence  of  the  bounded  model  uncertainty  ?(x,  t)  is 
compensated  and  the  asymptotic  stability  of  the  ASGSMO  is  guar¬ 
anteed.  As  the  switching  gains  are  adaptively  updated  in  response 
to  the  dynamics  of  system  state  errors,  the  ASGSMO  can  provide 
robust  tracking  capability  against  modelling  errors  and  significantly 
reduce  the  magnitude  of  chattering  in  the  SOC  estimation.  The  next 
section  will  present  the  experimental  results  to  validate  the  pro¬ 
posed  method  for  the  effectiveness  in  SOC  estimation. 

4.  Verification  of  ASGSMO  for  SOC  estimation 

The  configuration  of  the  ASGSMO  for  the  SOC  estimation  is 
illustrated  in  Fig.  14,  where  the  properly  defined  discharge  current 
profiles  are  simultaneously  applied  to  the  LiPB  integrated  with  Ah 
counting  module  and  the  ASGSMO  module.  The  LiPB  terminal 
voltage  is  sampled  and  fed  into  the  ASGSMO  to  generate  the  output 
tracking  error,  which  can  be  used  to  update  the  corresponding 
switching  gains  of  ASGSMO  so  as  to  compensate  the  modelling 
errors.  The  resultant  output  of  ASGSMO  module  is  the  estimated 
SOC,  which  is  concurrently  compared  to  the  true  SOC  directly 
generated  by  an  Ah  counting  module  to  show  the  accuracy  of  the 
SOC  estimation. 

The  feedback  gain  matrix,  I<  is  obtained  by  using  the  LQR 
method  with  Riccati  equation  as 

K=  [1.4244  1.4062  -0.1451  -0.2950]7"  (33) 

The  matched  model  uncertainties  are  determined  by  the  largest 
modelling  errors  after  comparing  experimental  LiPB  data  with 


X.  Chen  et  al.  /  Journal  of  Power  Sources  246  (2014)  667-678 


675 


Fig.  14.  Configuration  of  proposed  ASGSMO  system. 


model  data  so  that  the  bound  of  model  uncertainties  satisfy  |?(x,  t)\ 
<  0.2.  The  model  uncertainties  input  matrix,  r  is  defined  as 


r=  [1.4516  -0.0326  0.0122  0.0204]r  (34) 

And  according  to  an  algorithm  for  the  design  of  P  [32],  let 
Q.  =  2/4,  the  symmetric  positive  definite  matrix  P  of  the  LE  (22)  can 
be  determined  as 


P  = 


0.7075 

0.7032 

-0.0746 

-0.1500 


0.7032 

44.2880 

5.2992 

-6.2178 


-0.0746 

5.2992 

24.1883 

-0.6702 


-0.1500" 

17.6076 

-0.6702 

39.1701 


(35) 


The  positive  constant,  a  is  selected  as  0.5  to  satisfy  the  adapta¬ 
tion  speed  for  the  switching  gains  function.  To  validate  the  effec¬ 
tiveness  and  robustness  of  the  ASGSMO  for  the  SOC  estimation, 
three  types  of  discharge  current  profile  with  various  discharge  rates 
have  been  conducted  on  the  LiPB  at  the  approximately  room  tem¬ 
perature  of  25  °C.  As  the  battery  is  fully  charged,  the  initial  SOC  of 
the  LiPB  is  set  to  100%  before  loading  each  discharge  current  profile. 

In  order  to  demonstrate  the  better  performance  of  the  ASGSMO 
over  the  conventional  SMO  for  the  SOC  estimation,  three  current 
profiles  as  shown  in  Figs.  7-9  with  average  discharge  currents  of  1/ 


3Cn,  lCn  and  1.5Cn  have  been  used  to  conduct  experiment  on  the 
testing  platform. 

The  first  set  of  testing  data  is  obtained  by  using  the  CCD  profile 
with  discharge  rates  of  1/3 Cn,  lCn  and  1.5Cn,  namely  1.67  A,  5  A,  and 
7.5  A.  The  ordinary  SMO  and  the  proposed  ASGSMO  are  used  to 
estimate  at  any  initial  SOC  (e.g.  70%)  for  different  discharge  current 
rates.  As  an  example,  the  results  of  the  CCD  with  1  Cn  rate  are  shown 
in  Figs.  15  and  16,  respectively.  It  can  be  found  that  both  methods 
are  robust  to  incorrect  initial  SOCs.  With  relatively  large  constant 
switching  gains,  the  conventional  SMO  even  achieves  faster  SOC 
estimation  convergence,  but  the  significant  chattering  ripples  are 
observed  on  the  SOC  estimation  results  due  to  its  inability  of 
adaptively  updating  the  switching  gains.  Fig.  15  shows  its  estima¬ 
tion  errors  which  are  bounded  in  ±10%.  Comparing  with  the  SMO, 
the  ASGSMO  can  track  the  true  SOC  with  reasonable  convergence 
speed  and  small  chattering  ripples,  its  SOC  tracking  errors  are 
fluctuating  in  the  small  range  of  -2%  to  ±2%,  as  shown  in  Fig.  16. 
This  is  due  to  the  fact  that  the  switching  gains  of  the  ASGSMO  can 
be  adaptively  tuned  to  the  proper  low  levels  as  the  corresponding 
estimation  errors  decrease.  To  further  show  the  robustness  of  the 
proposed  ASGSMO  for  different  initial  SOCs,  Fig.  17  shows  the  RMSE 
of  the  SOC  estimation  errors  for  the  initial  SOCs  from  0%  to  100%  in 
the  step  of  10%,  they  are  all  within  3%. 


Fig.  15.  Ordinary  SMO-based  SOC  estimation  for  CCD  lCn. 


Fig.  16.  ASGSMO-based  SOC  estimation  errors  for  CCD  lCn. 
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The  second  and  third  sets  of  testing  data  are  obtained  from  city 
and  suburban  current  profiles  which  emulate  the  typical  driving 
cycles  as  a  driver  commuting  between  home  and  work  in  both  city 
and  suburban  settings  [31  ].  Both  city  and  suburban  current  profiles 
are  performed  with  average  discharging  rates  of  1/3  Cn,  lCn  and 
1.5Cn.  As  an  example,  the  SOC  estimation  results  of  two  observers  in 
city  current  profile  with  average  1  Cn  rate  are  shown  in  Figs.  18  and 
19,  respectively.  The  city  current  profile  as  given  in  Fig.  8  charac¬ 
terises  the  congested  traffic  conditions  in  city  with  repetitive  cycles 
of  acceleration  and  braking.  The  initial  SOC  is  also  set  to  70%.  For  the 
SMO-based  approach,  Fig.  18  illustrates  the  SOC  estimation  results 
and  the  estimation  errors  which  are  bounded  in  ±10%  with  sig¬ 
nificant  chattering  ripples.  For  the  ASGSMO-based  approach,  the 
SOC  estimation  errors  are  reduced  to  the  range  of  ±3%  and  -3%  for 
10-90%  of  SOC  with  minor  chattering  ripples,  as  shown  in  Fig.  19. 
Meanwhile,  the  ASGSMO  has  performed  superior  robustness  and 
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Fig.  19.  ASGSMO-based  SOC  estimation  errors  for  city  1  C„. 


capability  to  track  the  true  SOC  regardless  of  incorrect  initial  SOCs 
with  all  the  RMSEs  within  4%,  as  shown  in  Fig.  20.  As  another 
example,  the  SOC  estimation  results  of  both  SMO-based  and 
ASGSMO-based  methods  in  suburban  current  profile  with  average 
lCn  rate  are  shown  in  Figs.  21  and  22,  respectively.  This  profile 
characterises  the  smooth  driving  conditions  in  suburban  with 
minimal  stops.  With  the  same  initial  SOC  of  70%  as  the  city  profile, 
the  SMO-based  approach  gives  the  SOC  estimation  errors  within 
±10%  as  shown  in  Fig.  21  while  the  ASGSMO-based  approach  gives 
the  SOC  estimation  errors  in  the  range  of  ±3%  and  -3%  for  10-90% 
of  SOCs  with  nearly  smooth  SOC  estimation  results  as  shown  in 
Fig.  22.  The  RMSEs  of  the  SOC  estimation  for  different  initial  SOCs 
varying  from  0%  to  100%  in  the  step  of  10%  are  all  less  than  3.5%,  as 
shown  in  Fig.  23.  Table  3  summarises  the  RMSEs  of  the  SOC  esti¬ 
mation  for  each  type  of  current  profile  with  average  different 
discharge  rates.  It  can  be  seen  that  the  RMSEs  of  the  SOC  estimation 


SOC  [%] 


Fig.  18.  Ordinary  SMO-based  SOC  estimation  for  city  lCn. 


Fig.  20.  RMSEs  of  ASGSMO-based  SOC  estimation  at  different  initial  SOCs  for  city  1C„. 
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Fig.  21.  Ordinary  SMO-based  SOC  estimation  for  suburban  1C„. 


based  on  the  ASGSMO  is  always  lower  than  those  based  on  the 
conventional  SMO. 

To  further  illustrate  the  improved  accuracy  of  the  SOC  estima¬ 
tion,  the  probability  of  relative  SOC  estimation  errors  bounded  in 
5%  for  both  the  SMO  and  the  ASGSMO  under  the  above-mentioned 
three  different  current  profiles  has  been  evaluated  and  summarised 
in  Table  4.  It  can  be  seen  that  the  ASGSMO  has  performed  better 
SOC  estimation  with  approximately  99%  of  estimation  errors 
bounded  in  5%  than  the  SMO  with  only  approximately  83%. 
Therefore,  the  proposed  ASGSMO  can  provide  more  robust  and 
accurate  SOC  estimation  with  less  chattering  effects  than  the  con¬ 
ventional  SMO. 


SOC  [%] 


Fig.  23.  RMSEs  of  ASGSMO-based  SOC  estimation  at  different  initial  SOCs  for  suburban 
1C„. 


Table  3 

Comparison  of  SOC  estimation  RMSEs  based  on  SMO  and  ASGSMO  in  CCD/city/ 
suburban  profiles. 


Current  profile 

SMO  SOC 
estimation  RMSE 

AGSMO  SOC 

estimation  RMSE 

CCD  1.5Cn 

0.166 

0.026 

CCD  lCn 

0.154 

0.018 

CCD  1/3C„ 

0.092 

0.013 

City  1.5Cn 

0.247 

0.031 

City  1C„ 

0.192 

0.026 

City  1/3  Cn 

0.143 

0.017 

Suburban  1.5Cn 

0.233 

0.029 

Suburban  lCn 

0.184 

0.022 

Suburban  1/3  Cn 

0.132 

0.013 

5.  Conclusions 


The  adaptive  switching  gain  sliding  mode  observer  (ASGSMO) 
based  on  the  BECM  has  been  presented  to  estimate  battery  SOC  in 
EVs  and  it  is  able  to  compensate  modelling  errors  caused  by  the 
parameter  variations  of  the  BECM.  The  LiPB  is  utilised  to  conduct 
the  experiments.  The  parameters  of  the  BECM  are  extracted  from 
the  experimental  data  under  a  sequence  of  pulse  current  discharge. 
The  performance  and  effectiveness  of  the  SOC  estimation  based  on 
the  proposed  ASGSMO  have  been  verified  by  the  experimental  re¬ 
sults  of  a  LiPB  under  different  discharge  rates  of  various  current 


Table  4 

Comparison  of  probability  of  relative  SOC  estimation  errors  bounded  in  5%  based  on 
SMO  and  ASGSMO  in  CCD/city/suburban  profiles. 


Current  profile 

Probability  of  SMO 

SOC  estimation  error 
bounded  in  5%  (%) 

Probability  of  ASGSMO 
SOC  estimation  error 
bounded  in  5%  (%) 

CCD  1.5C„ 

79.36 

98.73 

CCD  lCn 

85.52 

98.95 

CCD  1/3C„ 

88.47 

99.21 

City  1.5C„ 

71.28 

98.76 

City  lCn 

76.34 

99.17 

City  1/3  Cn 

80.82 

99.32 

Suburban  1.5Cn 

77.76 

98.87 

Suburban  lCn 

81.92 

99.24 

Suburban  1/3  Cn 

84.29 

99.65 

Fig.  22.  ASGSMO-based  SOC  estimation  errors  for  suburban  lCn. 
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profiles.  It  shows  that  the  proposed  ASGSMO  has  a  robust  tracking 
capability  in  the  SOC  estimation  with  less  chattering  levels  and  high 
estimation  accuracy  in  the  comparison  of  the  conventional  SMO. 

Acknowledgment 

This  research  work  is  supported  by  Commonwealth  of  Australia, 
through  the  Cooperative  Research  Centre  for  Advanced  Automotive 
Technology  (AutoCRC),  under  the  project  of  Electric  Vehicle  Control 
Systems  and  Power  Management  (C2-801). 

References 

[1  ]  K.E.  Aifantis,  S.A.  Hackney,  R.V.  Kumar,  High  Energy  Density  Lithium  Batteries, 
Materials,  Engineering,  Applications,  Wiley-VCH  Verlag  GmbH  &  Co.  KGaA, 
Weinheim,  2010. 

[2]  T.B.  Reddy,  Linden’s  Handbook  of  Batteries,  fourth  ed.,  McGraw  Hall,  2010. 

[3]  S.  Piller,  M.  Perrin,  A.  Jossen,  Journal  of  Power  Sources  96  (1 )  (2001 )  1 13-120. 

[4]  V.  Pop,  H.J.  Bergveld,  P.H.L.  Notten,  P.P.L.  Regtien,  Journal  of  Measurement 
Science  and  Technology  16  (12)  (2005)  93-110. 

[5]  L.  Lu,  X.  Han,  J.  Li,  J.  Hua,  M.  Ouyang,  Journal  of  Power  Sources  226  (3)  (2013) 
272-288. 

[6]  K.S.  Ng,  C.-S.  Moo,  Y.-P.  Chen,  Y.-C.  Hsieh,  Journal  of  Applied  Energy  86  (9) 
(2009)  1506-1511. 

[7]  S.  Rodrigues,  N.  Munichandraiah,  A.K.  Shukla,  Journal  of  Power  Sources  87  (4) 
(2000)  12-20. 

[8]  H.  Blanke,  O.  Bohlen,  S.  Buller,  R.D.  Doncker,  B.  Fricke,  A.  Hammouche, 
D.  Linzen,  M.  Thele,  D.U.  Sauer,  Journal  of  Power  Sources  144  (6)  (2005) 
418-425. 

[9]  J.L.  Jespersen,  A.E.  Tonnesen,  K.  Norregaard,  L.  Overgaard,  F.  Elefsen,  World 
Electric  Vehicle  Journal  3  (5)  (2009)  2032-6653. 

[10]  C.C.  Chan,  E.W.C.  Lo,  W.X.  Shen,  Journal  of  Power  Sources  87  (11)  (2000) 
201-204. 

[11]  W.X.  Shen,  C.C.  Chan,  E.W.C.  Lo,  K.T.  Chau,  Journal  of  Energy  Conversion 
Management  43  (4)  (2002)  817-826. 

[12]  W.X.  Shen,  C.C.  Chan,  E.W.C.  Lo,  K.T.  Chau,  IEEE  Transactions  on  Industrial 
Electronics  49  (7)  (2002)  677-684. 

[13]  T.  Hansen,  C.-J.  Wang,  Journal  of  Power  Sources  141  (3)  (2005)  351-358. 

[14]  O.  Barbarisi,  F.  Vasca,  L.  Glielmo,  Journal  of  Control  Engineering  Practice  14  (3) 
(2006)  267-275. 

[15]  J.  Wang,  B.  Cao,  Q,  Chen,  F.  Wang,  Journal  of  Control  Engineering  Practice  15 
(5)  (2007)  1569-1576. 

[16]  A.  Vasebi,  M.  Partovibakhsh,  S.  Mohammed,  Journal  of  Power  Sources  17  (11) 
(2007)  30-40. 

[17]  B.S.  Bhangu,  P.  Benley,  D.A.  Stone,  C.M.  Bingham,  IEEE  Transactions  on 
Vehicular  Technology  54  (3)  (2005)  783-794. 

[18]  G.  Plett,  Journal  of  Power  Sources  161  (10)  (2006)  1356-1368. 

[19]  G.  Plett,  Journal  of  Power  Sources  161  (10)  (2006)  1369-1384. 

[20]  J.  Zhang,  C.  Xia,  Journal  of  Electrical  Power  &  Energy  Systems  33  (3)  (2011) 
472-476. 


[21  ]  W.  He,  N.  Williard,  C.  Chen,  M.  Pecht,  Journal  of  Microelectronics  Reliability  53 
(6)  (2013)  840-847. 

[22]  J.  Han,  D.  Kim,  M.  Sunwoo,  Journal  of  Power  Sources  188  (3)  (2009)  606-612. 

[23]  F.  Sun,  X.  Hu,  Y.  Zou,  S.  Li,  Energy  36  (5)  (2011)  3531-3540. 

[24]  F.  Zhang,  G.  Liu,  L.  Fang,  H.  Wang,  IEEE  Transactions  on  Industrial  Electronics 
59  (2) (2012)  1086-1095. 

[25]  J.  Yan,  G.  Xu,  H.  Qian,  Y.  Xu,  Energies  3  (10)  (2010)  1654-1672. 

[26]  I.  Kim,  Journal  of  Power  Sources  163  (10)  (2006)  584-590. 

[27]  X.P.  Chen,  W.X.  Shen,  Z.W.  Cao,  A.  Kapoor,  Australian  Journal  of  Electrical  & 
Electronics  Engineering  9  (3)  (2012)  225-234. 

[28]  M.  Chen,  A.  Gabriel,  M.  Rincon,  IEEE  Transactions  on  Energy  Conversion  21  (2) 
(2006)  504-511. 

[29]  H.  He,  R.  Xiong,  J.  Fan,  Energies  4  (3)  (2011)  582-598. 

[30]  X.  Hu,  S.  Li,  H.  Peng,  Journal  of  Power  Sources  198  (10)  (2012)  359-367. 

[31]  E.  Martinez-Rosas,  R.  Vasquez-Medrano,  A.  Flores-Tlacuahuac,  Journal  of 
Computers  and  Chemical  Engineering  35  (5)  (2011)  1937-1948. 

[32]  C.  Edwards,  S.  Spurgeon,  International  Journal  of  Control  59  (5)  (1994) 
1211-1229. 


Glossary 

A,  B,  C:  system  parameters  matrices 
Cn:  nominal  capacity  of  LiPB  (Ah) 

Cpe-  electrochemical  polarisation  capacitance  (F) 

Cpc:  concentration  polarisation  capacitance  (F) 

e(t):  state  estimation  error 

ey(t):  output  estimation  error 

I<:  system  feedback  gain  matrix 

Rin:  ohmic  resistance  (Q) 

Rpe-'  electrochemical  polarisation  resistance  (Q) 

Rpc ■  concentration  polarisation  resistance  (Q) 

Soc:  state  of  charge 

Soc:  estimated  state  of  charge 

Voc:  open  circuit  voltage  (V) 

Voc(Soc):  open  circuit  voltage  as  a  function  of  state  of  charge 
v  electrochemical  polarisation  voltage  (V) 

Vpe:  estimated  electrochemical  polarisation  voltage  (V) 

VjC:  concentration  polarisation  voltage  (V) 

Vpc:  estimated  concentration  polarisation  voltage  (V) 

Vt:  estimated  battery  terminal  voltage  (V) 

Vt:  battery  terminal  voltage  (V) 
x( t):  system  states 
y(t):  system  output  state 
a:  adaptive  law  parameter 
77:  Coulombic  efficiency 

AVoc,  A  Rin,  kRpe,  A  Cpe,  A  Rpc,  A  Cpc:  variations  of  circuit  parameters 
A/i,  A/2,  A/3,  A/4:  model  uncertainties 
T:  system  uncertainties  input  matrix 
h,  h,  h,  ^4-  eigenvalues  of  system  matrix,  A 
matched  uncertainty  function 
matched  uncertainty  bound 
v :  adaptive  switching  gains  function 


